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Abstract 

A simple and computationally efficient way of finding inverse avalanches for Abelian 
sandpiles, called the inverse particle addition operator, is presented. In addition, the 
method is shown to be optimal in the sense that it requires the minimum amount 
of computation among methods of the same kind. The method is also conceptually 
nice because avalanche and inverse avalanche are placed in the same footing. 
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The Abelian Sandpile Model (ASM), whose mathematical structure is first studied extensively 
by Dhar Q, is one of the few class of models of self-organized criticality in which a lot of interesting 
physical properties can be found analytically. The model consists of a finite number of sites labeled 
by an index set /. For each site i e /, we assign an integer hi called the local height to it. Whenever 
the local height of a site exceeds a threshold (which is fixed to for simplicity), the site is called 
unstable and it will transport some of its local heights (or sometimes called particles at that site) 
to the other sites in the next timestep by 

hj — > hj — Ay whenever hi > 0. (1) 

A is called the toppling matrix whose elements satisfies: 

An > Vie/, (2a) 



Ay < j, (2b) 

and 

]T Ay > Vie/. (2c) 
jei 

Toppling is repeated until all sites become stable again. The whole process of toppling is collectively 
known as an avalanche. The system is driven by adding unit amount of particles onto the sites 
randomly and uniformly after the system regains its stability. 

A system configuration, stable or not, can be regarded as a point in the space Z N where N is 
the total number of sites in the system. Both the addition of a particle and the toppling of particles 
in a site can be regarded as a translation in Z N PHI. The process of adding a particle to the site i 
together with the subsequent toppling it triggers (if any) can be viewed as a map between the set 
of all stable system configurations, and is denoted by aj UH]. 

Based on the observation that the final stable state is independent of the order of toppling in 
different sites, Dhar shows that a, o a, (a) = a.j o aj(a) for any stable system configuration a [|TJ. 
This is why we call the model 'Abelian". Using this commutative property, the total number of 
recurrence states in the model is shown to be det A ||, and we denote the set of all recurrence 
states by Q. 

Generalization of the ASM, known as the Generalized Abelian Model (GASM), has been made 
recently by Chau and Cheng H. In their model, the local heights hi and the elements in the toppling 
matrix A are real numbers instead of integers. Also, A satisfies only Eqs. (pa|) and (pbj). Arbitrary 
amount of particles are allowed to add to the system in possibly different locations all at the same 
time. Moreover, some special kind of configuration dependent triggering thresholds are used to 
determine the local stability of the pile (see ref. for details). In spite of the large differences 
between ASM and GASM, similar commutative properties between particle addition operators are 
found for recurrence system configurations. 

In both the ASM and the GASM, one can prove that for any pair of a e Q and a^, there exists 
a unique system configuration (3 e such that Bi(j3) = a 0J4]]. While the avalanche problem 
(i.e. the problem of finding a given (3 and a^) is straight forward and can be done very quickly in 
computer; the inverse avalanche problem (i.e. the problem of finding (3 given a and a^) is much 
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more difficult. (3 cannot be found, in general, by simply removing particle(s) from site(s) because 
particle removal operation can map a state out of the eventual phase space Q. Furthermore, the 
relationship a^ 1 = af et A_1 does not work very well for two reasons. First, det A is in general a huge 
number making the method computationally impractical. Second, in the event of GASM, det A 
may not be an integer and hence af etA is not well defined. 

The first computationally feasible method of finding inverse avalanches is proposed recently by 
Dhar and Manna using the so-called inverse avalanche operator by means of the burning algorithm 
||. However, the method works only on ASM with a symmetric toppling matrix. In this letter, we 
introduce a simple method to tackle the inverse avalanche problem for both ASM and GASM. To 
each particle addition operator, we find a nice and simple way to associate its inverse to another 
particle addition operator, called the inverse particle addition operator. Then we prove that the 
method is computationally optimal in the sense that it requires the least number of topplings among 
all the possible methods using the idea of inverse particle addition operators. 

For simplicity, we concentrate only on the case of Abelian Sandpiles in the discussions below. 
However, all arguments, after slight modifications, work equally well on GASM. We represent a 
system configuration a, stable or not, by a column vector of length N. In particular, the marginally 
stable state of the system (i.e. the one to which avalanche is triggered whenever particles are added 
to any one of the sites) is = (0, 0, . . . , 0). 

Given a particle addition operator a, we consider 

7 = (lj) je i = a (0) = (% ~ E fe n fc A fei ) ig/ G Q (3) 

where a« > is the number of particles added to site i, and n,i G Af for all i G I is the total 
number of toppling in site i triggered after the particles are added 0. For example, aj = Sij for 
the operation of adding a single particle to site i together with the subsequent toppling induced a^. 
The above definition works equally well when more than one particle is introduced to the system 
each time, and when they are introduced to different sites. 

Consider the operation of adding — ^ particles to site i for all z G / together with the subsequent 
toppling induced (if any). We denote this operation by a. This is a well defined operation sending 
system configurations from Q to Q because 7, < for all i. For any a = (ai) ig/ G fi, 

a(a) = (ojj - cij + T,k n k'A kj ) jeI G Q (4) 

for some G Z. Moreover, 

a o a(a) = (aj + T,k n k"^kj) jeI € ^ (5) 

for some nu" G Z. But by the remark between corollary 1 and corollary 2 in ref. M, we conclude 
that rik" = for all k G I and hence aoa(o) = a for all a G f2. Since the particle addition operators 
commute with each other, we can also conclude that a o a(a) = a. As a result, a = a^" 1 is the 
inverse particle addition operator corresponding to the particle addition operator a. 

We proceed to show that the above way of finding inverse avalanches is computationally optimal 
in the sense that the total number of toppling involve in the calculation is minimum among all the 
inverse particle addition operators (such as a detA_1 ). 

Suppose b' is another inverse particle addition operator corresponding to the particle addition 
operator a consisting of adding bi(> 0) particles to site % for allz G / together with the subsequent 
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toppling. Clearly, bj = J2k m k^kj — dj > for some e M @. Consider the addition of hi particles 
into site i l for all i at the same time to the system configuration 7. The resultant configuration is 
H = (J2k ( m k — n k) Afcj)j g /- Since ^(7) = a -1 (7) = 0, [i must either equal to or it is an unstable 
configuration which will eventually topple to 0. In either case, we can conclude that mi — rii > 
for alH £ / and the equality holds if and only if hi = — ji for all i £ I. 

For any a G Q, the introduction of 6j particles to site i for all i is equivalent to the addition 
of first — 7j particles to site i for all i and then Y^k{ m k — n k)^kt particles into site i for all i. So 
as compared to a, b' requires J2i{ m i — n i) more toppling to find the inverse avalanche. Thus a is 
computationally optimal. 

In summary, we have introduced a nice, simple and efficient way to find the inverse avalanche 
for both the Abelian Sandpile and the Generalized Abelian Sandpile by means of inverse particle 
addition operator. The method is conceptually nice because inverse particle addition operators are 
placed in the same footing as the particle addition operators. 
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